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1 Introdution 



A system of Schwingcr-Dyson equations (SDEs) is the infinite set of integral eqiia- 
tions containing, in principle, all the information about a model of quantum field 
theory. So far, there are no effective methods for the study of this infinite system 
as a whole, so it is necessary to truncate this system with a finite set of equa- 
tions. This truncation usually refers to either an expansion in a small parameter 
(examples of such expansions are the coupling-constant perturbation theory and 
the 1/A^-expansion), or a simulation (and rather, guessing) of some properties of 
the model (see, e.g., [1]). An initial truncation of the system determines further 
approximations of the complete system, which is thus attached to the leading ap- 
proximation. 

The system of SDEs is, in fact, the system of relations between derivatives of the 
generating functional of Green functions, resulting from the functional-differential 
SDE, which acts as a dynamical principle of the theory. If we approximate the gen- 
erating functional with the first n terms of an expansion in powers of a source, the 
system of SDEs can be approximated by a closed system of integral equations. This 
system defines the n-th term of a sequence of approximations, which forn — >■ oo ob- 
viously goes into the complete system of SDEs. A sequence of such approximations 
in the model of complex scalar field with interaction A((/)*(/))^ was considered in [2]. 
For the simplest nontrivial approximation (the "two-particle approximation") an 
asymptotic solution of the corresponding system in deep Euclidean region was ob- 
tained. In the strong-coupling region, this solution is free from Landau singularities. 
In [3], in the same approximation, the Yukawa model was considered, for which a 
self-consistent solution without Landau singularities in the Euclidean region also 
has been constructed. 

In this paper, in the same approximation, we consider a system of SDEs in the 
model of the complex scalar field (j> (phion) and real scalar field % (chion) with the 
interaction g<j)*(l)x in four dimensions. This model, also known as a scalar Yukawa 
model, is used in nuclear physics as a simplified version of the Yukawa model without 
spin degrees of freedom, as well as an effective model of the interaction of scalar 
quarks [4]. Despite its well-known imperfection associated with its instability [5] 
(or more precisely, the metastability [6] - [9]), this model, as the simplest model of 
the interaction of fields, often used as a prototype of more substantive theories to 
elaborate the various non-perturbative approaches in the quantum field theory. 

In the two-particle approximation, the system of SDEs for the scalar Yukawa 
model is a system of two nonlinear integral equations for propagators. The study of 
this system shows a change of the asymptotic behavior of propagators in the deep 
Euclidean region in a vicinity of a certain critical value of the coupling constant. 
For small values of the coupling the propagators behave as free, which is consistent 
with the wide-spread opinion about the dominance of perturbation theory for this 
super-renormalizable model. In the strong-coupling region, however, the asymptotic 
behavior changes dramatically - both propagators in the deep Euclidean region tend 
to some constant limits. 
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The existence of a critical coupling constant in the scalar Yukawa model was 
noticed by practically all authors who have investigated this model using diflFerent 
methods (sec, e.g., [10]- [16] and references therein). This critical constant is gener- 
ally regarded as a limit on the coupling constant for a self-consistent description of 
the model by some method. In our approach, however, the self-consistent solution 
for propagators exists also for the strong coupling, and the existence of the critical 
coupling looks more like as a phase transition in accordance with the general defi- 
nition of the phase transition as a sharp change of properties of the model with a 
smooth change of parameters (sec, e.g., [17]). 

The structure of the paper is as follows. In section 2, the general formalism of 
SDEs for the scalar Yukawa model is described and the necessary definitions and 
notations are given. A simplest non-perturbative expansion, namely the mean-field 
expansion is also considered in this section. The leading term of this expansion cor- 
responds to the chain summation for a chion propagator. The problem of restoring 
of crossing properties for this expansion is discussed, and the existence of a critical 
value g1 of the coupling is shown. For g"^ > get the asymptotics of inverse chion 
propagator in deep Euclidean region becomes negative, which leads to the appear- 
ance of a Landau singularity and the associated violation of the self-consistency of 
the method. 

In section 3, a construction of n-particle approximations for the system of SDEs 
is given, and the two-particle approximation is considered in detail. The system of 
SDEs in the two-particle approximation is reduced to a system of nonlinear integral 
equations for propagators. The study of the system at large Euclidean momenta 
shows that for the equal masses of fields a self-consistent solution exists for each 
value of the coupling, and there is a critical coupling value, which separates the weak 
and strong coupling regions with a different asymptotic behavior. For the unequal 
masses, the weak and strong coupling regions are separated by the segment of the 
intermediate coupling without self-consistent solutions. 

In section 4, a three-particle approximation for the model is briefly considered, 
and the crossing-symmetry problem for two-particle amplitude is discussed. A 
discussion of the results is contained in section 5. In the Appendix, the exact 
solution of linearized integral equation in the strong-coupling region is obtained for 
the case of equal masses. 

2 SDEs and the mean-field expansion in the bilo- 
cal source formaUsm 

2.1 Preliminaries 

We consider the model of interaction of a complex scalar field (j) (phion) and a real 
scalar field x (chion) with the Lagrangian 

C = -d^<P*d^<p - mlr<l> - lid^xf - yx' + fl-^Vx (1) 
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in a four-dimensional Euclidean space {x S -E4). The coupling constant g has a 
dimension of mass . 

The generating functional of Green functions (vacuum averages) is the func- 
tional integral 

G{v,j) = J D{(j),(j)*,x)expl^j dx/:{x)-j dxdy^*{y)ri{y,x)^{x)+j dxj{x)x{x)^. 

(2) 

Here, 77 is a bilocal source of the phion field, ^ j is a single source of the chion field. 

The translational invariance of the functional integration measure leads to the 
functional-differential SDEs for the generating functional G. In terms of the loga- 
rithm Z = log G these equations are 



5^Z 



-+ 



6Z 6Z 



5r]{y, x)5j{x) 6j{x) Sr]{y, x) 



SZ 



Sr]{y,x) 



+ j dxi r]{x, 



Xi) 



6Z 



6z_ _ r 



dxi Dc{x - xi) - g dxi Dc{x - xi) 



I' 



dZ 



6r]{xi,xi)' 



_5Z 

Here Dc = {n"^ - d"^)-^ . 

Equation (4) allows us to express all Green functions with chion legs in terms 
of functions that contain phions only. Thus, the differentiation of (4) over 77 gives 
us the three-point function 

g / dzi Dc{z - zi) Z2 



+d{x-y), 
(3) 
(4) 



V{x,y\z) 
where 



dj{z)6r]{y,x) 



r)=j=0 

X y \ _ 
x' y' 



Zl Zi 

x y 



(5) 



(6) 



v=j=o 



S'q{y',x')Sr]{y,x) 

is the two-particle phion function. The differentiation of (4) over j with taking into 
account equation (5) gives us the chion propagator: 

S^Z f 

Dc{x-y)+g^ / dxidyiDc{x-xi)Z2 



D{x~y) 



^j{y)Sj{x) 



Xi Xi 

yi yi 



Dciyi-y), 
(7) 



etc. 



Excluding with the help of the SDE (4), a differentiation over j in SDE (3), we 
obtain at j = the SDE for the generating functional: 

S'^Z SZ SZ 



g'^ J dxi Dc{x - 



Xl) 

SZ 



Sri{xi,xi)Srj{y,x) Sri{xi,xi) 5ri{y,x) 
SZ 



Sr]{y,x) 



dyiv{x,yi) 



Sv{y,yi) 



+ S{x -y)=0, 



(8) 



The formalism of the bilocal source was first elaborated in the quantum field theory 
by Dahmen and Jona-Lasinio [18]. We consider this using as a convenient choice of the 
functional variable. 
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which only contains the derivatives over the bilocal source rj. 

Sequential differentiations of this equation give us the infinite system of SDEs 
for Green functions. For our purposes, it will need the first three equations of this 
system. Switching off the source in equation (8), we obtain the equation 



(m^ - dl)A{x -y) = 5{x - y) + j dxt D^x - xi) 



X y 

Xi X\ 



(9) 



Here = mo - ^A(a; = 0) and 



A(a; -y) = - 



5Z 



5r]{y,x) 



(10) 



7;=0 



is the phion propagator. 

A differentiation over ry gives us the second equation: 



{m''-dl)Z^(^l, y,^-g^ j dx,D,{x-x,)Z^ 



Xi Xi 

x' y' 



A{x -y) + 



^ I Xi Xi \ 

dxiD^{x-xi)Zs\ X y ] = 6{x - y')A{x' - y). {11) 
J \x' y' j 

The repeated differentiation over r/ gives one more equation: 

I X y \ ^ / xi xi \ 

{m^-dDZsl x' y' \-g^ dx^D,{x-x^)zA x' y' \ A{x - y) 
\ x" y" I J \ x" y" j 



~^9^ j dxi Dc{x — xi) Z4 
+g'^ j dxi Dc{x - xi)Z2 
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Xl 


Xl 


\ 
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x' 


y' 
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y" 


J 



Xl Xl 

x' y' 



x" y" 



Here ^„ = gf 



X y 

is the n-particle phion function. 



-y) + 

^{x-y")Z.[^l„ 

x' ^x\y' ^y"\.(yi) 



77=0 



2.2 Kernel and inverse kernel 

Integral equations with the kernel 

X y 



Kab 



x' y' 



5{x-x')5{y-y')-g'^5{x' -y') j dxiAa{x-xi) D^{xi-y')Ab{xi-y) 

(13) 
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will be repeatedly considered below . Here Aa and A(, are the given functions. It 
is easy to verify that the inverse kernel has the form 



Kb ( ^/ ^=6{x-x')6{y-y')+6{x'-y') J dx,Aa{x-x,) f{x,-y')A^ 



b{xi-y), 
(14) 



J dx,dy,K-J ) K^, ^} ^) ) = 5{x - x')S{y - y') (15) 



I.e. 



(no summation over a and 6!). 
In equation (14), 

f-\x-y) = ^D-\x-y)-K,{x-y), (16) 

where 

Kb{x-y) = Aa{x-y)At,{y-x). (17) 
Note also the useful property of the inverse kernel: 



5^ j dxi D^{x - xi) KJ p 



f{x-y')S{x'-y'). (18) 



2.3 Mean-field expansion 

The mean- field expansion for the generating functional (see [2], [3] and references 
therein) 

^^^(o) + ^(i)+... (19) 
is based on the leading approximation: 

/ dxiDc{x-xi) — — r+(m'-a^)— -+ / dyir]{x,yi) — r+S{x-y) = 0. 

J dr](xi,xi) dT]{y,x) Or](y,x) J or}{y,yi) 

(20) 

The NLO equation is 

g / dxi Dc[x - xi) — — +5 / dxi Dc[x - xi) 



Sr){xi,xi) Sr]{y,x) J ' 5ri{xi,xi) Sr]{y,x) 

- + y -^^^ ^(^'^^^ ^^K^ \ ^-^-^ - ^^) ^,(.„.,)^,(,,x) ' (^^^ 

etc. 

Equation (20), after switching off the source, gives us the leading-order phion 
propagator. In the momentum space 

Ao'(p) = m2+p2. (22) 
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A differentiation of the leading-order equation over 77 gives us tlie equation for 

(0 

2 



the leading-order two-particle function z!f '' , which can be written as: 



I dx^dy,Koo ( ) (5 p ) = - y')Mx' - y), (23) 

where Kqq is kernel (13) at = A;, = Aq. In correspondence with equation (14), 
the solution of this equation is 

y,^=Aoix-y')Ao{x'-y) + 

+ j dxidx2Ao{x - xi)Ao{x' - X2)fo{xi - a;2)Ao(a;i - y)Ao{x2 - y'), (24) 

where f^^ = g~^D~^ — Lq, and Lo = AqAq is the scalar loop. Equation (7) after 
simple calculations with taking into account above formulae gives us the leading- 
order chion propagator Dq. In the momentum space 

D^\p) = l,^+p^-g^Lo{p^). (25) 

The repeated differentiation of the leading-order equation over r/ gives us the 
equation for the leading-order three-particle function whose solution is 

X y \ f / \ I x\ y\ 

Zf\ x' y' \ = \ dx^dy^K^n I I ) Z30 x' y' \ . (26) 



x" y" ) ■> \xx yx ; , ^„ y„ 



Here 



V/<i»iO.(i-m)4°'( 5 P)4°'(^^, », ) }-|a,'«i",s'«/|,(27) 

Note the following property of the three-particle function resulting from equa- 
tion (18): 



/. I Xi Xi 

/ dxiDc{x- xi)Z'f^ X y 
J \x' y' 



1 Xi 


Xl \ 






y 


\-l 




y' ^ 





Xl Xl 

dxi fo{x - Xl) Z30 { X y \ ■ (28) 
x' y' 

The other functions of the leading approximation can be calculated in the same 
manner. 
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Note that in contrast to the functional derivatives of Z over single source j, the 
higher derivatives of Z over bilocal source are not connected parts of correspond- 
ing many-particle functions. Thiis, the two-particle phion function Zi is not the 
connected part Z^ of two-particle and related to it by the formula: 



(29) 



A characteristic feature of many- particle functions of the leading approximation 
is their incomplete structure in terms of crossing symmetry. The Bose symmetry 
of the theory dictates the crossing symmetry of the connected part Z^: 

'.)-<:i »')• « 

It is easy to see that zf^'^ (this is the second term on the rhs of equation (24)) sat- 
isfies the first equality in (30), but breaks the other two. This apparent discrepancy 
is a feature of many non-perturbative approximations. It is inherent, for example, 
to the Bethe-Salpeter equation in the ladder approximation. We face similar prob- 
lems in the two-particle approximation, which will be considered below. To restore 
the missing crossing symmetry of the leading approximation, it is necessary to look 
at the next order. 

Calculations in the next order quite similar to the above. Equation (21) with 
the source being switched off gives us the equation for the first correction Ai to the 
propagator. The differentiation of equation (21) over the source gives the equation 
for the NLO two-particle function z'^^ , which can be written as an equation with 
kernel Kqq, and taking into account the above formulae the solution is 

^2'^ ( i'' ) " ( 5 ^' )^/ '^^i'^^2Ao(x-xi)/o(a;i-a;2)Ao(a;i-y)Z; 



21 

(31) 



X2 X2 

x' y' 



where 



X y 
"''^ X' y' 




y 
y' 



= -9^ j dxidyiAo{x - yi)Dc{yi - xx)z[ 

+Ao(a; - 2/')Ai(a;' -y)- ^A,{x = 0) j dx^Aoix - a;i)zf 

J dxidx2Ao{x - xi)De{xi - X2)Ai{xi - y)Z^^^ -'^^ ^ . (32) 

From these equations, together with equation (28), it follows that Zj^^ contains the 
term 

j dxidx2Ao{x - a;i)Ao(a;' - a;2)/o(a;i - a;2)Ao(a;2 - y)Ao(a;i - y') = 
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"j=<"(: i), (33) 

which restores the missing crossing symmetry of the LO two-particle function. Such 

restoration of the crossing symmetry is typical for non-pcrturbativc expansions in 
the formalism of bilocal source (for similar examples in other models see [19], [20]). 

The above equations contain divergent integrals and require a renormalization. 
According to the standard recipe wc introduce the renormalized propagators 
and Dj. and impose on them the normalization condition^ 



A-i(O) 
and 



ml, 



dA- 



dD- 



= 1 (34) 

p=0 



dp^ 



= 1. (35) 

p=0 

In the leading approximation of the mean-field expansion the phion propagator is 
given by equation (22), which implies that it is sufficient to replace Aq — >■ A^., ml — >■ 
m^, i.e., the renormalized leading-order phion propagator is 

A;\p)=ml+p^. (36) 

For the renormalized chion propagator according with equation (25), we obtain in 
the leading approximation: 

D;\p')=l^l+p^-g^Lr{p^), (37) 

where L,(p2) = Loip'') - Lo{0) - p''L'o{0). 

At p"^ -)> 00, 

D-\p') = (1 - ^;;^)p' + 0(log 4)- (38) 

As can be seen from this equation, the asymptotic behavior of the phion propagator 
in the deep Euclidean region is self-consistent in the region of weak coupling which 
in this case is determined by the condition < = 967r^m^. When g^ > g^, 
the asymptotics of the inverse propagator becomes negative, which leads to the 
appearance of a Landau singularity for the chion propagator in the Euclidean region, 
which in turn means a violation of basic physical principles. As was mentioned in 
the Introduction, the existence of similar restrictions on the value of the coupling 
constant in the model has been noted by many authors. Some authors (see, e.g.. 



^For the easement of the following calculations we choose the normalization point at 
zero momenta. Such a choice is not a big deal for simple calculations of the mcan-ficid 
expansion, but it is highly significant for more complicated calculations in a more pithy 
two-particle approximation (see below). 
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[14], [15]) believe self-evident that the presence of this kind of critical constants 
reflects the metastability of the model. It seems to us, however, that the presence 
of such a singularity means firstly the non-applicability of the calculation method 
in the strong coupling region and the need for more meaningful non-perturbative 
approximations. 



3 Two-particle approximation 

The system of SDEs, generated by functional-differential equation (8), is an inflnite 

set of integral equations for n-particlc phion functions Zn = 6"Z/Sri'"\ji=o. The first 
three SDEs are equations (9), (11) and (12). The nth SDE is the (n— l)th derivative 
of SDE (8) with the source being switched off and includes a set of functions from 
onc-particlc function (phion propagator) to the (n + l)-particlc phion function. In 
order to obtain a sequence of closed systems of equations, we proceed as follows. We 
call "the n-particle approximation of the system of SDEs" the system of n SDEs, in 
which the first n — 1 equations arc exact and the nth SDE is truncated by omitting 
the (n-h l)-particle function. It is evident that the sequence of such approximations 
goes to the exact set of SDEs at n — >■ oo. The one-particle approximation is simply 
equation (9) without Z2. This approximation has a trivial solution which is a 
free propagator. The two-particle approximation is a system of equation (9) and 
equation (11) without Z^: 

{m'-dl)Z2 ( ^, y, yg'l dx,D,{x-x,)Z2 ( J ^ ) A{x-y) = S{x~y')A{x 

(39) 

This system of equations will be the object of the present investigation. 

An alternative method to obtain this system is a modification of the mean-field 
expansion of section 2. In this modified mean- field expansion, equations (9) and 
(39) are the basic equations of the expansion. A construction of such modified 
mean-field expansion was performed for the model \{(j)*4>)'^ in work [2] and easily 
can be extended to the considered model. 

In framework of the two-particle approximation. Z2 can be expressed as a 
functional of A. In fact, equation (39) can be written as 



X 


y 




xi yi \ 








y' 



dx,dy,Ki I ]ZA f; \=A,{x-y')l\{x'-y), (40) 



where K = Kab at = Ac = (m^ — d"^) ^ and A;, = A (see (13)). In correspon- 
dence with equation (14), we obtain 

Z^l^l, y,^=A,{x-y')A{x'-y) + 

+ j dxidx2 Ac{x - xi)A{x' - X2) f{xi - a;2)A(a;i - y)Ac{x2 - y'), (41) 
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where / is given by equation (16) and 

L{x-y)^A,{x-y)A{y-x). (42) 

Equations (7) and (9) with the taking into account above formulae and equation 
(18) give us the system of equations for the chion and phion propagators. In the 
momentum space, this system has the form 

D-l(p2)=^2+p2_^2^(p2) 

where 

Hp') = / + 9)A(<Z), K{p^) = I ^ A,(p - q)D{q). (44) 

The renormaUzation of equations of the two-particle approximation is similar 
to the renormalization of equations of the mean-field expansion with the imposition 

of normalization conditions (34) and (35) on the renormalizcd propagators. The 
system of renormalized equations for the chion and phion propagators is: 

\ A-\p^) =ml+p-'- g'Kr{p% 

where 

Lr = L{p^) - L(0) - p^i'(O), Kr = K{p'^) - K{0) - p'^K'{0). (46) 

The system of equations (45) is a system of non-linear integral equations, which 
is quite difficult for the analytical study. In the study of the asymptotic behavior in 
the deep Euclidean region, one can make the approximation, which greatly simplifies 
calculations, namely to replace Ac by the asymptotics at p"^ — > oo, i.e. by the 
massless propagator 1/p^. This massless integration approximation is quite usual 
in investigations in the deep Euclidean region [21]. 

Using the well-known formula 



d'^q $(g2) _ 1 



(27r)4 {p-qY 16^2 



4/ ^q^)q^dq^+ ^q^) dq' 
P Jo Jp^ 



(47) 



gives us the system of equations for propagators in the massless integration approx- 
imation: 

' A-Hp^) =m^ + il- ^)p^ + ^ if dq^D{q-){l - ^) 
D~\p-)=^,- + {l-^)p- + ^^fdq^A{q-){l-^^). 

Here and below we omit the index r, bearing in mind that all quantities are renor- 
malized. 

Note that the system of equations (48) is symmetric with respect to the change 

A^ D, w? ^ fj,^. (49) 
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With taking into account this symmetry, we introduce the dimensionless variables: 



(Here m = y/m?, = \fi^-) 

In terms of these variables, the system of equations (48) takes the form 

ru(t) = (^-A)t+i + 2A/J4i;^(i-^) 



(50) 



The normalization conditions (34) and (35) for the dimensionless functions u and 
V have the form 

u(0) = u(0) = 1; u{0) = ^, viO) = -. (51) 

m fi 

From the system of integral equations (50) one can easily go to the system of 
nonlinear differential equations for the functions u and v: 

l^M = 2(^-A) + f. 

The boundary conditions for this system are the normalization condition (51). 

If V (or u) increases in the absolute value at f — )• oo, then for A ^ m/ /x, A 7^ f^/m, 
we obtain from the system (52) the asymptotic behavior in the deep Euclidean 
region: 

r« = (^-A)i+o(iogt) 

\v = {f-X)t + 0{logt). ^^'^> 

This asymptotic behavior is self-consistent in the weak-coupling region at 
A < min{^, ^} ( < min{/x^, m^}) and corresponds to the asymptotically- 
free behavior of propagators: 

A ~ 1/p^ D - 1//. (54) 

Outside the weak-coupling region one should distinguish three cases: 
^2 = ^2 ^ ^2 g^^^ ^2 ^ j^2 

At /i^ = m^, the system of equations (50) takes the form: 



v = {l-X)t + l + 2Xj^^{l-'^). 



(55) 



An iterative solution of this system gives us u^") = v^""^ for any finite number n 
of iterations; hence, u = v (provided the convergence of the iteration procedure). 
Unfortunately, it is difficult to prove the convergence of this nonlinear iteration 
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scheme. Nevertheless, we wih assume the equahty of u and v in the study of the 
case /U^ = rm?, leaving the possibility of an alternative for a subsequent study. 
At u = the system of equations (55) is reduced to the single integral equation 

.= (1-A)t + 1 + 2A^^|;^(1-^), (56) 
which, in turn, can be reduced to the differential equation 

^(..) = 2(1-A) + -. (57) 

In the strong-coupling region A > 1, equation (57) has the positive exact solu- 
tion ^ 

Uexact = Ua = ^ _ ■ (58) 

Note that Uq is a solution of the integral equation 

uo = {l-X)t+- ^+2a/— (1-i), (59) 

1 - ^ Jo Uo t 

which differs from equation (56) by the inhomogeneous term, and this difference is 
small for large t in comparision with the leading part (1 — \)t of the inhomogeneous 
term. This points that uq is an asymptotic solution of equation (56) at t ^ oo. In 
fact, performing the linearization: 

1 1 Ul 

U = Uo + U\, -~ 5-, (60) 

U Uo Uq 

we obtain for ui the linear differential equation 

_(iui) = -a^ui, (61) 

where 

a = (A -1)^1. (62) 
The solution of the linearized equation (61) is 

Ul = (vli Ji (2a Vt) + ^2 n (2a Vt)) . (63) 

Here Ji Yi are the cylindrical functions.^ Therefore, ui = o{uq) at large t, and 
constant uq is really the asymptotic solution in the strong-coupling region. 



*For exact solution of the linearized integral equation see Appendix. 
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At the critical value A = 1, the differential equation (57) has the exact solution 



Uexact = U= X/—. (64) 

u is the solution of integral equation 

/■* dt' t' 

u = 2 (65) 
Jo u t 

which differs from equation (56) at A = 1 by the inhomogeneous term. It is an 
indication that u is the asymptotic solution of equation (56) in the critical point 
A = 1, since this inhomogeneous term is o{u) at t — )• oo. Indeed, performing the 
linearization 

_ 1 1 Ml 

u = u + ui, - ^ - - (66) 
we obtain the linearized equation for u\: 

^(i«i) = -|«. (67) 
which is an Euler equation, and the real solution of this equation is 

If. . logi , logt\ 
Ui = —^[Ax sin — ^ + A2 cos — ^ . (68) 

Therefore, ui = o{u) at t — >■ 00. 

Thus, the asymptotic behavior of propagators at the critical point A = 1 has 
the form 

A - l/v/p2, D ^ 1/^. (69) 

For unequal masses, we will consider taking into account the symmetry (49) only 
the case fj? > m?. In this case, there are two critical couplings: Ad = m//U and 
Ac2 — l^/'rn, which determine the weak-coupling region (A < m/ii), the intermediate 
region (m//z < A < /u/m) and the strong-coupling region (A > /u/m). 

At A = Aci = m/n, as in the weak-coupling region, there is a self-consistent 
solution with the asymptotic behavior for large t of the form 

u ~ ^ 1, V ^ ^ ^ log t. 70 

yum /Lt^ — 

In the strong-coupling region A > Ac2, the system of differential equations (52) 
has the positive exact solution 

Xfi Am 

Uexact = Uo = , Vexact = Vq = . (71) 

AjjL — m Am — jjL 
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{uq, vq} is the solution of the system of integral equations 



|.o = (^-A)i + ^+2A/3'f(l-4) 



(72) 



As in the above case of equal masses, this system differs from the system (50) by 
inhomogeneous terms, and this difference is 0(1) at large t, while the inhomoge- 
neous term is 0{t). This fact indicates that {uq, vq} is an asymptotic solution at 
t — > 00. 

The linearization of the system of differential equations (52), completely analo- 
gous to the linearization of equation (57), leads to linear equations, the solution of 
which, as well as for the case of fj,^ = m^, is expressed through cylindrical functions 
and has the order of o(l) for large t, which confirms the conclusion that {uq, vq} 
is the leading asymptotic solution. 

In the intermediate region m//i < X < fi/m, there are no self-consistent solu- 
tions. 

In fact, if either u, or ?; increase in absolute value, then the system of differential 
equations (52) dictates asymptotic behavior (53), i.e. the function v is negative at 
large t and, therefore, the chion propagator has a Landau-type singularity. 

If the function v is positive and bounded, i.e. < w < Vmax, then from the 
integral equation for u (see (50)), we obtain u > (p/m — A + X/vmax)t + 1, i.e. the 
function u increases, and this increase again leads to the infinite increase (in the 
absolute value) of v in contradiction with the assumption. 

A similar analysis shows that at A = Ac2 = /u/m a self-consistent solution also 
does no exist. 

Thus, in the case of unequal masses the critical point becomes a segment of 
critical values Ad < A < Ac2, in which self-consistent solutions are not available. 

4 Three-particle approximation 

The two-particle approximation as well as the leading term of the mean-field ex- 
pansion has an incomplete crossing structure of two-particle function. For the two- 
particle approximation the violation is even more significant because it affects the 
disconnected part. However, in the same way as for the mean-field expansion, this 
problem is solved by considering the next approximation. The next three-particle 
approximation is described by a system of three equations, the first two of which 
are equations (9) and (11), and the third is equation (12) without Z^. 

Just as in the two-particle approximation, Z2 can be expressed as a functional 
of A, in the three-particle approximation Z3 can be expressed as a functional of Z2 
and A: 




(73) 
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where 



X y 
x' y' 
x" y" 



A,(x- - y")Z2 



X y 
x" y 



+ 



+0^ j dxidx2 Ac{x - xi) Dc{xi - X2)Z2 





X2 




' x" 


y" \ 




y' . 




. ^1 


y J 



-l^x' ^x",y' ^y"y (74) 



and is defined by equation (14) at = A^ = A. Using this expression 
and the properties of the inverse kernel K~^, we obtain an equation for Z2 in the 
three-particle approximation, which can be written as 



X y 



x' y' 



X y 
x' y' 



dxidyi Ac{x-xi)f{xi-yi)A{xi-y)Z2 
wher Z2 is the functional of Z2 and A: 

Z2 \ I I 

\x y 



yi yi 

x' y' 
(75) 



y ^ -.A,{x-y')A{x' -y) + 

+ j dXidX2 A,{X - X,)f{xi - X2)I^A,{X2 - y')Z2 2 ^ + Ae(X2 - y)Z2 ^ + 



+5^ j dyidy2Ac{x2 - y2)Dc{y2 - yi 



Z2 
+Z2 



yi yi 

xi y 

yi yi 

x' y' 



Z2 

Z2 



x' y' 

y2 X2 

xi y 

2/2 X2 



+ 



(76) 



If we approximate in Z2[Z2, A] the dependence on Z2 by the two- particle approxi- 
mation, i.e. we consider the iterative solution, in which Z2[Z2,A] Z2[Z|^'*, A], 
then we can show, using equation (9), that in this approximation the correct cross- 
symmetric structure is restored for the disconnected part and also for the connected 
part of the two-particle function. 



5 Conclusions 

The main result of the work is the finding of the change of asymptotic behavior in 

the scalar Yukawa model in the framework of two-particle approximation. The sys- 
tem of SDEs in the two-particle approximation has self-consistent solutions in the 
weak-coupling region (where a dominance of the perturbation theory is obvious) 
and in the strong-coupling region. The field propagators in the strong-coupling 
region asymptotically approach constants. It is not something unexpected, if we 
look at the results of studying the strong-coupling region in models of quantum 



16 



field theory. In particular, the well-known result in this direction is the conception 
of the ultra-local approximation (also known as the "static ultra-local approxima- 
tion"), considered in the paper CaianicUo and Scarpetta [22]. This exactly soluble 
approximation is based on removing of kinetic terms 5^ in the Lagrangian. As a 
result, all the Green functions are combinations of delta functions in the coordinate 
space that arc constants in the momentum space. Of course, this approximation is 
very difficult for a physical interpretation. Nevertheless, it can be considered as a 
starting point for an expansion in inverse powers of the coupling constant, i.e. as a 
leading approximation of the strong-coupling expansion 

In this connection, it is noteworthy that the obtained solutions of the two- 
particle approximationin in the strong-coupling region for large Euclidean momenta 
tend to constants. Such behavior asymptotically corresponds to the ultra-local 
approximation. On the other hand, our solutions are free from the interpretation 
problems, since for the small momenta they have the quite traditional pole behavior. 
This indicates that this approximation seems to be adequately described as the 
weak-coupling as the strong-coupling too. 

In the case of equal masses, of fields a self-consistent solution of SDEs in two- 
particle approximation exists for any value of the coupling, including the critical 
value. At = g\, the asymptotics of propagators are i.e. the asymptotic 
behavior is a medium among the free behavior at 5^ < g\ and the constant- 
type behavior in strong coupling region > g^- A sharp change of asymptotic 
behavior in the vicinity of the critical value is a behavior that is characteristic 
for a phase transition. This phase transition is not associated with a symmetry 
breaking, and in this sense is similar to the phase transition of "gas liquid". The 
weak-coupling region can be roughly classified as the gaseous phase and the strong- 
coupling region, where a kind of localization of correlators exists - to the liquid. 
This analogy, of course, is a quite shallow. A type and characteristics of this phase 
transition can be defined as the result of a detailed study using methods of the 
theory of critical behavior. 

In the case of unequal masses, two critical values of coupling exist. In the 
interval between them, it is not possible to construct a solution with self-consistent 
ultraviolet behavior. Reasons for the existence of such an interval is currently 
unclear. The existence of such intermediate values of coupling can somehow reflect 
the metastability of the model. On the other hand, it may be an artifact of the 
two-particle approximation. In any case, a detailed study of the three-particle 
approximation of Section 4 will help to shed a light on this issue. 

■'The ultra- local approximation and the strong-coupling expansion based on it are dis- 
cussed in detail in the book of Rivers [23]. For later works, see [24], [25]. For the ultra- local 
approximation in the bilocal source formalism, see [26]. 
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Appendix 



A linearization of nonlinear differential equations in the strong-coupling region has 
a sense, in general, only for large t. However, in the particiilar ease of eqiiation (57) 
at A > 1, this linciarization available throughoiit the range of the t, starting from 
zero, i.e. a linearization of the integral equation (56) is possible. We can use this 
fact to estimate the values of constants Ai and A2 in equation (63). Performing the 
linearization of integral equation (56) by formula (60), we obtain for ui the linear 
integral equation 



1 f t' 

ui = -- — I- -a' / dt'ui{t'){l- -), 
^ Jo * 



(77) 



where a is defined by equation (62). Equation (63) tells us to deal with the Ansatz 
in the form 

ui = ^Ji{2aVi), (78) 
yt 

i.e. A2 = 0. The value of constant Ai is determined by the direct substitution 
of the expression (78) into integral equation (77). Calculating integrals with the 
known formulae (see, e.g. [27]), we obtain 

Thus, the solution of the linearized approximation is 

u = uo + ui = -^-- ^— ^Ji(2av^). (80) 

A — i (A — l)avt 
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